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Abstract 

We have investigated the denaturation bubble life time for a homogeneous as well as for a 
heterogeneous DNA within a Poland- Scheraga model. It is shown that at criticality the bubble life 
time for a homogeneous DNA is finite provided that the loop entropic exponent c > 2 and has a 
scaling dependence on DNA length for c < 2. Heterogeneity in the thermodynamical limit makes 
the bubble life time infinite for any entropic exponent. 
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Introduction. DNA denaturation has a long history . With variation of temperature 
or the solvent pH value double-stranded DNA denatures, yielding regions of single stranded 
DNA, until the double strand is fully molten at a melting temperature T^. The melting 
temperature is between 60° — 100"C depending on the ratio of AT and GC pairs in DNA. 
However, a local base opening is possible even at room temperatures \^. This fact explains 
the biological relevance of local denaturation. In particular, it is relevant to the transcription 
process where a DNA is partially opened for RNA polymerase to read the template. A DNA 
bubble consists of two flexible single strands with fluctuating size. The bubble size varies 
from a few broken pairs well below T^, up to several hundreds closer to T^. A well-known 
model of DNA melting is Poland-Scheraga model 

In the Poland-Scheraga model, the partition function of a single bubble has the form: 



Z{y) = aoexp 
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where y is the bubble size, i.e the number of broken AT or GC base pairs, 6{y,T) is a free 
energy released by a broken base pair, (3 = 1/T and ctq = 6"^°*^^ is the statistical weight for 
initiation of a loop. In physiological conditions it is very small ~ 10^^ which means that 
the formation of a loop is a rare effect and one can consider a single loop instead of many 
loops. The function e{y,T) describes the heterogeneity of the loop. 

new experimental 



possibilities opened by molecule fluorescence correlation spectroscopy(FCS) [l3l]- It enables 
one to monitor the bubble dynamics in a real time. One of important problem in the bubble 
dynamics is the length dependence of the averaged bubble life time. Previously it was 
reported [11] that for c > 1 the bubble life time is finite and in the thermodynamical limit it 
does not depend on the length both at critical and non-critical temperatures. In the present 
article we study this problem within the Poland-Scheraga model in a systematic way. We 



get results different from those obtained in 11|. 



Initial Relations. We assume that e{y, T) has the form 

eiy,T)=eiT)+r]iy,T) (2) 

where s{T) ~ Tc — T is an average value and 7]{y,T) is a random function describing the 
randomness of DNA sequence consisting of AT and GC pairs. The factor (1 + takes into 
account the lack of entropy of the loop compared to a single strand with the same length. 



The free energy of the bubble correspondingly takes the form 



ry 

F{y) =^0+ e{y)dy + cTln{l + y) 
Jo 
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The relaxation time of a bubble is much longer (10 — lOOyUs) than the relaxation times of 
polymeric degrees of freedom. Therefore it is reasonable to assume that the bubble dynamics 
takes place on the background of thermodynamical equilibrium for polymeric degrees of 
freedom and is governed by a stochastic Langevin equation 

dy r^dFiy) 
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where q{t) is a thermal gaussian white noise 



<q{t)q{t')>=2DT8{t-t'), <^(t)>=0 
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and D is the diffusion coefficient in the base pairs unit. 

Using Eqs.([3}lMj can obtain the Fokker-Planck equation for a probability density 
function (PDF) (jj 
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It is convenient to introduce new units HI] and rewrite Eq.(j6])in a simpler form 
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where we redefine e s/2T, t — » 2DTt, v = c/2 and ri{y,T) — > ri{y,T)/2T. Eq.(I71) is a 
generalization of the corresponding Fokker-Planck equation for a homogeneous DNA 
At the denaturation transition point the number of opened base pairs is large y ^ 1. Then 
if one ne glec ts the entropic term in Eq.([7]) i^/{y '^Xh problem is reduced to the Sinai 
diffusion [17] of a particle in random environments fiol. [isl. However we will keep it because 
as we will see below it plays important role in the denaturation bubble dynamics. 

Homogeneous DNA. In order to simplify the Fokker-Planck equation Eq.(I7]) we introduce 
a new function P 

P{y,t) = {l + y)-'P{y,t). (8) 

Substituting Eq.([8]) into Eq.([7]), one obtains an equation for P 

dP{y,t) d[e + r]{y)]P{y,t) 



dt dy 
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As it follows from Eq.([9]), far distant from the melting temperature e S> ?7, therefore 
heterogeneity will not play an essential role in bubble dynamics. We will consider the 
bubble dynamics at criticality e = where the role of heterogeneity is more important. 
First we consider a homogeneous DNA case. In this case ri{y) = and the Fokker-Planck 
equation Eq.([n]) is reduced to a Scrodinger equation in imaginary time 

dPiy,t) 



dt 

with Hamiltonian 



HP{y,t) (10) 



Taking into account the initial conditions and using Eqs.fl8l)- f|TTl) one can represent the PDF 
of bubble in the following form 

i^j $^e-^"*^„(?/)^.(yo). (12) 

^ ' n=l 

Here we assume that at t = the bubble size is ?/o,and where En and \l/„ are the eigenvalues 
and eigenfunctions of Hamiltonian Eq. ffTTl) . Note that the PDF obeys the following initial 
condition 

P(2/, 0,2/0,0) = 5(2/ -2/o)- (13) 

following from the completeness of eigenfunctions 

For the finite DNA the eigenfunctions should satisfy the boundary conditions \l/ri(0) = 
\1/„(L) = 0. Note that Eq. (fT2l) is also correct for heterogeneous DNA, because the Fokker- 

n 

Planck equation can be reduced to a Scrodinger equation in general case and the corre- 
sponding potential energy term will be random. At denaturation temperature the number 
of opened base pairs is large ?/ ^ 1, and one can neglect the centrifugal term in Eq. flllj) . 



Note that using the potential l/y^ instead of l/(y + 1)^ as in ll|] is incorrect because l/y^ 
contains an artificial divergence at ?/ = 0, which affects the long time asymptotics. Instead 
we can use the wave functions of quantum particle in a box with infinite walls at and L 
with eigenvalues En = ■n'^v? /2LF' 

^n(l/) = /|sin^. (14) 

As follows from the form of En and Eq. (fT2|) . for t ^ L^, one can be restricted by the first 
term in the sum and PDF will be exponentially small . For t <^ L^, the sum in Eq. (|T2l) can 



be substituted by an integral which is taken exactly 
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For 1/ < 1 it is safe to neglect the unity in entropic factor(see below). Using Eq. (fT5l) one 
can calculate the DNA bubble survival probability in the thermodynamical limit L —>■ oo, 
C{yo, t) = /o°° dyP{y, yo, t) 

r(2 - u) 



C{yo,t) 



r(V; 
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(16) 



which is correct provided that u < 2 and y^ <^t <^ Li^. Scaling law in Eg. (1161) differs from 



the law obtained in 
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ll| . C{t) ~ t . Bessel functions instead of sin were used 







when obtaining this law. The problem with this Bessel functions is that they do not satisfy 
the correct boundary conditions at L and are unstable relative to perturbations. 

Let us estimate the relative contributions of kinetic and potential energy terms of Hamil- 
tonian Eq. ffTTj) . Using Eq. ffT^ one can estimate the contribution of kinetic energy term 
as ~ 71^ / LF' and potential energy term at denaturation as ~ 1/L^. The typical value of 

~ L^/t follows from the cut off in Eq.( fT2|) due to the exponential function. Therefore 
the kinetic energy term contribution is of order ~ 1/t which, for the times t -C L^, is much 
larger than the potential energy contribution. Henceforth, the neglecting of potential energy 
term is justified. 

One can find the average bubble life time using the first passage time distribution function 
W{t) [jj 



tW{t)dt, 



(17) 



where the first passage time distribution function is expressed by the survival probability 
W{t) = —dC{t)/dt. Decomposing the integral (fTTIl into integrals over [0, L^] and [L^,oo] 
one obtains r ~ L}^'^ . This scaling is a transitional intermediate behavior between finite 
closing time at T < and finite opening time dXT > T^.. The bubble life time is maximal 
at criticality. For u > 1 the bubble life time is finite even at criticality. 

In the existence of two absorbing boundaries the first passage time describes the time at 

n 

which the brownian particle reaches one of the boundaries for the first time [15]. Applying 
this to our problem means that r is either opening or closing bubble time at the criticality. 
So, in our consideration the critical value of loop entropic ex pon ent is z/ = l(c = 2) in 



consistent with the statics of denaturation transition 



19 
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2l| . Remind that for c > 



2 the denaturation transition is of first order which imphes a finite bubble hfe time in 
dynamics. For c < 2 the denaturation transition is of second order and the bubble life time 



imit. 
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ll| introduces another critical value 



correspondingly is infinite in the thermodynamical 

Note that the scaling result C{t) ~ t^^^'^^'^ from 
for dynamics u = l/2(c = 1) which does not coincide with its value in statics. 

Let us explain why for investigation of asymptotic behavior of C{yo, t) the set of functions 
Eq. ffT^ is preferable compared to the set $fc(y) = {kyY^^ Jii2+v{ky), where J are Bessel 
functions which are exact eigenfunctions of the Hamiltonian Eq. (fTT|) with potential energy 
term v{v + l)/2y^ and were used in 
Hamiltonian Eq. fllip has the form 



,lll| . Exact solution of the Scrodinger equation with the 



= vVTT [Ci4+i/2(zA;(y + 1)) + c^K.+y^my + 1))] , (18) 

where k'^ = 2E, I and K are modified Bessel functions, Ci, C2 are arbitrary constants which 
should be found from the boundary conditions \E'fc(0) = "^kiL) = 0. Consider the asymptotic 
behavior of Eq. ffTSj) . For small ky <^ 1, taking into account that \1'(0) = and expanding 
Eq. lITS!) on y, one obtains "^kiy) ~ ky. Such a behavior differs from {kyY^" behavior in 
$fc(?/) case. The extra {kyY dependence is caused by the artificial divergence of potential 
l/y"^ at y = 0. The correct potential is smooth at y = 0. Note that small ky behavior of 
wave functions infiuences the long time asymptotics of survival probability C{t). For large 
ky ^ 1, using asymptotic expressions for / and K one can find "^kiy) ~ sin ky. Therefore 
the box functions ([HD, in contrast to ^k{y), correctly approximate the exact solution Eq. (fT8|) 
both for large and small y. 

Heterogeneous DNA. At critical temperature the number of opened base pairs is large 
y ^ 1. Therefore one can neglect terms proportional to ri{y) /{y + 1) and l/(y + 1}^ in Eq.Q 



and reduce the problem to the diffusion of particle in random environments [18[. Taking 
into account Eq.([H]) and using the method of images to satisfy boundary condition at y = 0, 
for the averaged over the disorder PDF of heterogeneous bubble, in the thermodynamical 
limit L — > 00, one has 

Piy, t, yo) = f Psiy, t, yo), (19) 

a 



Psiy,t,yo) 



In't 



InH 
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(20) 
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and 

Q(X) = -y izir e-^''^l(2'"+^)^ (21) 

m=0 

Here a is the variance of the random field ri{y) 

<r]{y)r]{y' >=2a6{y-y'). (22) 

Our approach differs from the one of [lOj since we include entropic effects which are 
neglected there. Using Eqs. fll9l) - fl2Tl) one can calculate the survival probability of the bubble 
C{yo,t) = P{y,t,yo)dy. After integration on y, we get 

^ ^^0, t) ^3 2^^2m + l (5(2m + l)2)i-- 



X 



^2 



e 



^g5{2m+l)2p^^ - z/, 5(2m + 1)2) - e-'^(2,n+l)2p^^ _ (23) 

where 'y{x,y),r{x,y) are incomplete gamma functions, T{x) is the proper gamma 

2 

function, (5 = "^^^2° and we assume that u < 1 . Let us now consider the asymptotic be- 
havior of C{yo, t). Substituting asymptotic expressions for 7, T into Eq. (|T6l) . one obtains for 
(5> 1. 

09 °° ( 1 \m 

c{y^^t)-4ll7pkv.^ (24) 

TT-^ ^-^ [2m + I r 

m=l ^ ' 

In the large time limit 5 <^ 1, expanding Eg. (1161) over 5 we obtain 

m=l 

As it follows from Eq.( l25|) the sum diverges for z/ > 0. However this divergence is artificial 
and caused by the expansion on b. Indeed, the sum (1231) converges for large m. Therefore 
the divergent sum in Eg. (1251) should be cut on some value M, at which the expression under 
the initial, correct sum fl23|) becomes of order of unity ~ 1. Hence the divergent sum 
(!25|) can be estimated as follows 

i6r(i - v)b^^-M''- i6r(i - v)b^^-i^ 

C{yo,t} . (26) 

Note that the peculiarity in prescaling factors in Eqs. (fT6!) and (!26l) for u > 1 is associated 
with the fact that we neglect the unity in the entropic factor (1 + y)"*^ and introduce an 
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artificial divergency at y = which turns out to be important for u > 1. Actually the 
scaling laws Eqs. f|T6|) and fl26|) are caused by the behavior at large y therefore with some 
non-peculiar prefactors they will be correct also for u > 1. Using Eq. fl26p one can find the 



first passage time 



15| distribution function for large times 



dC{t) 1 , ^. 



It follows from Eq. (1271) that the integral determining the bubble life time in denaturation 
r ~ tW{t)dt diverges. It should be cut at some large time scale. As was mentioned 
above, the representation (IT^ is correct also for a heterogeneous DNA. In this case En and 
\E'„(?/) are the eigenvalues and eigenfunctions of a random Hamiltonian consisting of terms 
vri{y,t) / {y + 1) and ri'{y,T)/2. Taking into account that ri[y,T) and ri'{y,T) are limited 
when y ^ oo and treating them perturbatively with the set Eq. (fT4l) one can show that, 
similar to the homogeneous case En ~ n^/L"^ when n — > 0. Using the same reasoning as in 
the homogeneous case one finds that PDF of heterogeneous DNA is also exponentially small 
for t ^ Cutting off the integral on the scale one finds the dependence of the bubble 
life time on DNA length in denaturation 

L2 



T ~ 



(28) 



(/nL)3+2^/3 

It follows from Eq. (!28|) that in the thermodynamical limit L — oo the bubble life time 



22| that the 



is infinite for any u. This result correlates with the exact result in statics 
heterogeneity makes the order of denaturation transition higher than one. 

In summary, we have considered the bubble life time problem in denaturation for homo- 
geneous as well as for heterogeneous DNA. For a homogeneous finite length DNA a scaling 
behavior is found when the entropic index c < 2 . For c > 2 the bubble life time is finite 
in the thermodynamical limit. Heterogeneity strengthens the L-dependence and leads to a 
significant increase of the bubble life time. It becomes infinite in the thermodynamical limit 
for any value of entropic index c. At non- critical temperatures a DNA bubble in a finite 
time is closing when T < and is opening when T > T^. 12|. Our results are consistent 
with the statical picture of denaturation transition. 
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